Roughness of a Tilted Anharmonic String at Depinning 
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We consider the discretized model of a driven string with an anharmonic elastic energy, in a two 
dimensional random potential, as introduced by Rosso and Krauth 0. Using finite size scaling, we 
numerically compute the roughness of the string in a uniform applied force at the critical depinning 
threshold. By considering a string with a net average tilt, we demonstrate that the anharmonic 
elastic energy crosses the model over to the quenched KPZ universality class, in agreement with 
recent theoretical predictions. 

PACS numbers: 74.60.Ge, 02.70.Tt, 05.60.Cd 
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Recently, Rosso and Krauth (RK) reported new 
simulations of roughening at the depinning threshold of 
a driven one dimensional string in a two dimensional (2d) 
random potential. Introducing higher order anharmonic 
terms to a quadratic elastic energy for the string, RK 
found a value for the roughness exponent £ ~ 0.63, in 
contrast to the values £ <~ 1.2 found in earlier simula- 
tions 2], and recently theoretically 0], using a purely 
quadratic energy. RK noted that the value £ ~ 0.63 had 
previously been found in some cellular automata models 
for depinning [4j. 

In a subsequent work 5] RK (with Hartmann) noted 
that that when an average tilt is applied to the string, the 
anharmonic terms break the rotational invariance present 
in the quadratic model, thus suggesting that the anhar- 
monic terms might cross the model into the quenched 
KPZ universality class, previously introduced by Kardar 
|3| to explain the anisotropic depinning observed in the 
automata models. Simulations |jj of a continuum model 
with the quenched KPZ term found C — 0.61 ± 0.06, 
consistent with the automata models and with the an- 
harmonic model of RK. Most recently, a functional renor- 
malization group calculation by Le Doussal and Wiese Q 
argued that the quenched KPZ term can indeed be gener- 
ated, not only by the anisotropic disorder considered by 
Kardar, but also by the anharmonic elastic energy terms 
introduced by RK. 

A key prediction of Kardar's for anisotropic depinning 
is that the roughness exponent for a tilted interface will 
differ from that of an untilted one; for a tilted string 
in 2d, he predicted the exact value of (tut = 1/2. In 
this report, by computing the Ctiit of RK's model for the 
first time, we offer a direct numerical demonstration that 
their model does indeed belong in the quenched KPZ 
universality class. 

Our model is the same as that of RK. We take for the 
energy of the string, 
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{V(i,hi)- fhi + E^hi+x-hi)} , (1) 



0, ... L on a discretized lattice, V(i,j) is an uncorrelated 
random Gaussian potential with zero average and unit 
variance, / is a uniform external driving force, and E c \ is 
the elastic energy of deforming the string. V(i,j) is taken 
periodic on an L x L system size. In their work, RK used 
periodic boundary conditions. Here, to model a tilted 
interface with net slope s, we use boundary conditions 
hr, = ho + sL. Defining the height relative to a uniformly 
tilted line, Shi = hi — si, so that Sh^ = 8ho, we can 
rewrite Eq.Q) in terms of the Shi and recover the same 
model as RK except that the elastic term now has the 
form, 
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We now carry out simulations of Eq. Q , using the elas- 
tic term of Eq.(J2Jl. Using the same algorithm [9j as RK, 
we consider slopes s = and s = 1 for one of the specific 
cases studied in 0, 



E el (A) = A 4 /16 



(3) 



We compute the interface roughness W for a system of 
length L, 



L-1 



(4) 



where hi, is the integer height of the string at position i 



where Shi ^ s the relative height at site i of the critical 
string at depinning, Sh c is the average relative height of 
the critical string, [. . .] represents an average over many 
realizations of the random potential V(i,j), and C is the 
roughness exponent. We also compute the disorder aver- 
age of the critical force, f c . 

Our results for string roughness W 2 vs. L, averaged 
over 500 disorder samples (for L = 2048 we use only 200 
samples), are plotted in Fig.^ For s = 0, our numerical 
values agree with those in Ref.Q. The straight lines on 
the log-log plot indicate the power law relation, W 2 ~ 
L 2< » , and the difference in slopes indicate clearly different 
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FIG. 1: Roughness W 2 vs. system size L, for strings of net 
slope s = and s = 1. The lines are the best fits to W 2 ~ L 2 ^ 
using system sizes L = 128 — 2048. 
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value of 1/2 predicted by Kardar, and thus verifies that 
the anharmonic model of RK is in the quenched KPZ 
universality class. 

Finally, in Fig.|31we plot the critical force f c as a func- 
tion of system size L for the tilted and untilted strings. 
We see clearly that the critical forces approach different 
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FIG. 3: Critical force f c vs. system size L, for strings of net 
slope s = and s = 1. The lines are guides to the eye. 



values as L increases, another signature of anisotropic 
depinning. 
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FIG. 2: Roughness exponents £ for the tilted (s = 1) and 
untilted (s = 0) strings, as obtained by fitting the data of 
Fig. □ to W 2 ~ L 2C using data for sizes L min to 20484. 

roughness exponents for the tilted (s = 1) and untilted 
(s = 0) strings. 

To determine the values of the exponent £, we fit the 
data in Fig.nto W 2 ~ L 2c =, using system sizes from L m i n 
to L = 2048. We plot the resulting values of £ vs. L m i n in 
Fig-Hi f° r ^min = 4 to 256. We see that as L m - m increases, 
the values of £ decrease and saturate to a constant value 
characterizing the roughness in the asymptotic large L 
limit. Using the results from fitting with L m i n = 128 we 
find for s — the value £ ~ 0.63 ± 0.01. This value is 
used to plot the solid straight line in Fig. ^ an d agrees 
with the value found by RK. For s = 1, however, we 
find the value ( = 0.52 ± 0.01. We use this value to plot 
the dashed line in Fig. ^ Given that ( for s — 1 still 
shows a small systematic decreases as L m in increases, we 
believe our value is in excellent agreement with the exact 
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